Optoacoustic solitons in Bragg gratings 
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Optical gap solitons, which exist due to a balance of nonlinearity and dispersion due to a Bragg 
grating, can couple to acoustic waves through electrostriction. This gives rise to a new species of 
"gap-acoustic" solitons (GASs), for which we find exact analytic solutions. The GAS consists of an 
optical pulse similar to the optical gap soliton, dressed by an accompanying phonon pulse. Close to 
the speed of sound, the phonon component is large. In subsonic (supersonic) solitons, the phonon 
pulse is a positive (negative) density variation. Coupling to the acoustic field damps the solitons' 
oscillatory instability, and gives rise to a distinct instability for supersonic solitons, which may make 
the GAS decelerate and change direction, ultimately making the soliton subsonic. 

^ ' PACS numbers: 42.81. Dp, 42.70.Qs, 43.25.+y, 72.50.+b 
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' Introduction. — Light and sound tend to have widely disparate frequencies, wavelengths, and velocities. As a conse- 
quence, interaction between optical and acoustic signals is often weak. The interaction is important, e.g., for Brillouin 
scattering, diffraction of light by ultrasonic waves Brillouin scattering in optical fibers — coupling of longitudinal or 
transverse acoustic modes to optical waves through electrostriction — has been analyzed extensively, though generally 
not for solitons 0, [^, 0, S 0] ■ Sound waves acting on trains of optical solitons in fibers have been studied 0, 13 , 
but these interactions were essentially between acoustic continuous waves and their optical counterparts, subject to 

■ periodic modulation. Brillouin scattering has been employed in acousto-optic filters, in which light diffracts off an 
acoustic grating 0. Brillouin scattering in fibers can lead to pulse compression Q or formation of trains of very 
narrow optical solitons [ld| . Interaction of phonons with optical nonlinear Schrodinger solitons, with nonlinearity due 

\ to electromagnetically- induced transparency, was analyzed in Ref. Reference [l^ considered solitons in Bragg 

■ gratings, with electrostriction the only nonlinearity; the correct soliton solution is a limiting case of the results below. 
\ In this work, we consider light and acoustic waves in a cubic (x^^'') nonlinear medium with a Bragg grating, a fixed 

OO ■ spatially periodic variation in the index of refraction [l^ . Bragg gratings are typically written in the cladding of an 
\ optical fiber, though they have been induced by copropagating continuous waves Fiber Bragg gratings without 

■ electrostriction are known to support op tical gap solitons, nonlinear solitary waves with frequencies in the band gap 
CD ' created by the Bragg grating [l3, 1^, Gap solitons can have arbitrarily small velocity, which has motivated great 



interest in them for stopping light, though creating very slow gap solitons experimentally has been difficult. Early 
T\ \ experimental realizations of gap solitons had velocities approximately half the speed of light [lij . Various methods 
^ ■ have been proposed for creating slow gap solitons: building up slow gap solitons directly by Raman amplification 
[13] , obtaining quiescent solitons from collisions of fast-moving gap solitons [l8| , and retarding fast solitons by passing 
them through gradually varying ("apodized") Bragg gratings [19[. In a recent experimental breakthrough, optical 
\ gap solitons were slowed to a sixth the speed of light using apodized Bragg gratings [l^l • 

Electrostriction, compression of the medium due to variations of the intensity of light, allows light to drive acoustic 
waves; conversely, acoustic waves feed back into the optical field through the dependence of the refractive index 
on the material density. Gap solitons can be arbitrarily slow; at velocities close to the speed of sound, an opto- 
acoustic resonance may dramatically enhance the interplay between light and sound. We derive a generalization of 
the standard optical gap soliton, including an acoustic component. The dynamics of these "gap-acoustic solitons" 
(GASs) arc significantly different from optical gap solitons without the acoustic dressing when the velocities are slow, 
and slightly different when velocities arc large. Our results demonstrate that if an optical gap soliton is slowed to a 
velocity on the order of the speed of sound, the inherent clectrostrictive effects can do the work of slowing the soliton 
to significantly below the speed of sound, even creating stopped light. The results also suggest that GASs may be 
controlled by continuous acoustic waves or pulses. 

We assume one transverse electromagnetic mode and one acoustic mode, either because the waveguide has exactly 
one transverse optical mode and one acoustical mode, or because coupling to other modes is negligible. We have 
derived effectively l-|-l-dimensional coupled mode equations using standard methods 13, 21| for light in a waveguide. 
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adding a phonon field and acousto-optic coupling induced by electrostriction [l|, 0, H, 0i Hi 0, B : 

= ifcottt + iMz + + '^'^^u^i^'^ (Xsl^P + XxK'P) M + XesWM, (la) 

= ik'^vt-iv, + K*u+ ^""^"^o/^) (xx|Mp + XskP)t^ + XcsW^', (lb) 
= u^tt-r^t,, + A(|M|2 + b|')... (Ic) 

Here z is the coordinate along the waveguide and t is time; u(z,t) and ^(z, t) are slowly varying envelopes of the 
counterpropagating electromagnetic waves {total electric field E{z^ t) = u{z, t) exp[i(fco2: — oJat)] + w(z, t) cxp[— i(A:o-z + 
uj^t)] + c.c.}, with carrier frequency ljq and wave numbers rkkiujo) = ±n(ujQ)ujo / c, n{Lo) being the index of refraction, 
fcp = dk/duj\i^=^g the reciprocal group velocity, and k the Bragg reflectivity. The self- and cross-phase modulation 
coefficients arc Xs and Xx, and A is the effective cross-section area of the waveguide modes. The acoustic field wlz, t) 
is the density of the medium, offset by a constant Wq. Equations ^ assume that the acoustic wavelengths arc much 
larger than the optical one, so sound waves are not subject to Bragg reflection. The form of the electrostrictive 
interaction involves the average square electric fleld {E'^) ~ 2(|mP -|- jwp) (see Eqs. (9.3.11-12) in Ref. [H). The 
electrostrictive coefficient A, which comes from the dependence of the energy density on the index of refraction 
n(a>. Wo + w), is A = (2t:)-^ {Wo + w)ndn/dw 0; in a waveguide, the coefficient may be strongly affected by the 
modal structure 0,13, Si- The other electrostrictive coefficient, Xca = (luq / c) dn/ dw , comes from local changes in the 
wave vector[/c = n{u!o, Wq + w)lu/c], which depends on the index of refraction, which in turn depends on the acoustic 
field [H, 0, H [ill- The speed is sound is denoted by I3s. Phonon viscosity T is included for completeness, though 
current Bragg waveguides are not long enough for it to be significant. Mass 



M ^ A I w{z, t)dz 

and number of photons 



are conserved; momentum 



A 



N=^-^AI (M^ + lv^dz 



dz . 



where r(z, t) = w{z\ t)dz' is an acoustic potential, is invariant for F = 0, but decays when viscosity is non-zero. 
Equations ([T|) have resemblances to the Zakharov system 2^, 23 1: both contain dispersive high-frequency waves. 



coupled to nondispersive acoustic modes. The dissimilarity is the source of dispersion — a Bragg grating rather than 
intrinsic material dispersion. Looking at Eqs. |T]) as a Zakharov system with a Bragg grating helps to identify 
realizations in additional physical systems: surface and bulk acoustic waves in solids, where the grating (for the 
surface mode) is a surface superstructure; coupled vibron and acoustic waves in polymer molecules, with the grating 
(for the vibrons) provided by an intramolecular structure; and coupled Langmuir and ion-acoustic waves in plasmas 



(the original physical Zakharov system), with a grating in the form of a dust crystal 24 1 



Soliton solutions. — We found a family of GAS solutions to Eqs. ([T]) with zero phonon viscosity (F = 0), 



~ i/k7(1 + /3 fcg) a sinQ sech(C sin Q — z(5/2) exp[i ^(C) — i T cos Q] , (2a) 



'k*7(1 - Pk'g)a sinQ sech(CsinQ -|- iQ/2) exp[i 6'(C) - ircosQ] , (2b) 
Al^a^l 47sin^ Q 

(zc) 



/32 - /32 cosh(2C sin Q) + cos Q ' 
with T = ^\K\[t/k'o~ {(3 k'o)z],C = l\^\{z-f3t), J =[l-{pk'oy']-^^^ and 
Qf.. 4|a|2(/3fc^) /27r(wo/c)^ , XcsA \ _i 

^(C) = "1 /^,,N2 1—A + 3^ tan '[tanh(Csmg)tan(Q/2)] , (3a) 



l-(/3fc^,)2 V koA (3^-/3^ 
_ 2^(a-o/c)2 / l + WKT] , 2xc.A 



koA V ^'l-{l3k',rj (/3?-/32)[l-(/3A:^) 



FIG. 1: (Color online) Three components of a supersonic gap-acoustic soliton (right- and left-traveling envelopes of the elec- 
tromagnetic waves and the acoustic field), for k'o = k = 27r(wo/c)^(fco^)~^Xx ~ Xcb = 1, Xs = Xx/2, (3 = 0.25, /3s — 0.2. F = 0, 
and A — 0.1. Dashed, dotted, and solid lines show, respectively, the real and imaginary parts and absolute values of the fields. 

The solitons are characterized by two parameters: f3, the soliton velocity, and Q, which takes values < Q < tt and 
determines the soliton's full width at half-maximum. (|k|7 sin Q)"^ cosh^^(2 + cosQ). peak intensity, and frequency 
[in the rest frame, 7|K|(fcg)~^ cos Q]. These solutions also hold for quiescent solitons (/3 = 0) with non-zero phonon 
viscosity (T > 0). The velocity may take any value up to the group velocity of light, < 1/fco, except for a 
velocity-spectrum gap at and above the speed of sound, ^ [/?s,/3ci], with 

2{k',)^x.-Xs 2 V V2(fc^)'Xx-Xs 2) 2^(c.o/c)2 Xx - Xs ' ^' 

For the typical case, Xx = 2xs > 0, the negative sign applies. For subsonic {\f3\ < f3s) solitons [supersonic solitons 
{fis < /?cr < \P\ < 1/^0)] ' density perturbation around the soliton is positive [negative], i.e., compression 
[rarefaction]. The light components {u,v) vanish as the velocity approaches the speed of sound from below, \P\ — > /3s; 
both light and sound (w) amplitudes diverge when the velocity approaches /3cr from above. Figurc[l]shows a moderately 
supersonic soliton, with frequency in the middle of the bandgap {Q ~ tt/2), and velocity /3 = (1.25)/3s = 0.25. 

With electrostriction, A 7^ 0, Eqs. ([T]) do not admit solitons without an acoustic component, but purely acoustic 
waves (solutions to the free D'Alembert equation for w) are possible. Equations ([T]) do not admit stable breathers, 
as oscillations of a localized mode will generate emission of acoustic waves, dissipating the oscillations. 

Stability. — We tested soliton stability by numerical simulation of Eqs. ([1]), using a split-step fast Fourier transform 
scheme [21|. Simulations were carried out systematically for three values of the soliton coefRcientjQ = 7r/3 (in the 
middle of the top half of the band gap), which, for gap solitons without electrostriction (A = 0) [l^, is well inside 
the stable region; Q = 7r/2 (mid-point of the band gap), which is stable but close to the instability border; and 
Q = 27r/3 (in the middle of the bottom half of the band gap), which has an oscillatory instability [2^. We took ten 
values of the electrostrictive coefficient A, ranging over four orders of magnitude, and the limit with zero instantaneous 
Kerr nonlinearity Xs = Xx = 0. The self- to cross- phase ratio was Xs/Xx = 1/2. The speed of sound was taken to 
be Ps = 0.2/fcg, which is faster than physically realistic, but more clearly illustrates the dynamics. Initial soliton 
velocities were taken over a range from zero to twice the speed of sound. To impose consistent perturbations, all 
numerical simulations had initial light amplitudes 1% greater than those of the exact soliton solutions. 

Like gap solitons without electrostriction [2^, GASs are subject to oscillatory instabilities, which grow until they 
destroy the solitons. However, electrostriction decreases the growth rate of the instability, as shown in Fig.[2l Similarly, 
numerical simulations showed that the growth rate of the oscillatory instability decreases for velocities closer to the 
speed of sound. A general trend is that larger phonon components produce greater damping of the oscillatory 
instability. The damping may be due to emission of acoustic waves, and that the GAS's acoustic component creates 
an effective trapping potential which suppresses disintegration of the optical components. 

All solitons that are stable without electrostriction remained stable with clcctrostiction A 7^ 0, if the velocity was 
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FIG. 2: (Color online) Time evolution of the peak power of the GAS, [maxz(|it(z, J)p + for electrostrictive coefflcients 

A = 0.0001, 0.001, 0.01, and 0.1. The solitons are taken with zero velocity and intrinsic parameter Q = 2 7r/3, which is unstable 
in the standard gap soliton model (A = 0). The speed of sound is Ps ~ 0.2. 





FIG. 3: (Color online) Instability of a supersonic GAS, for Q — tt/3 and initial velocity /3 — 0.25, which is 125% the speed of 
sound /3a ~ 0.2. The electrostrictive coefficient is A — 0.001. 



subsonic < f3s- AH supersonic (/3s < /?„■ < \P\ < l/^o) GASs were found to be unstable if numerically integrated 

for long enough. This supersonic instability is qualitatively different from the oscillatory instability, and is unknown 

for optical gap solitons without electrostriction. The supersonic instability may be connected to the fact that for 

slightly supersonic solitons, the soliton momentum decreases with velocity above the speed of sound 4^ < 0, which 

I — I _ .... 

is a known instability criterion [23[. Supersonically unstable GASs tend to retain their integrity through growth of 

the instability. The closer is the supersonic soliton's velocity to the speed of sound, the sooner the instability takes 

effect; this may be explained by the divergence of the momentum slope ^ at the critical velocity. The changes in 

velocity arc fairly abrupt, and arc accompanied by emission of phonons. which carry off momentum. The solitons may 

change speed and direction a few times before eventually settling to a stable subsonic form. When the electrostrictive 

coefficient is larger, and the phonon field larger, the supersonic instability tends to be stronger: it occurs sooner, and 

the soliton is more likely to be destroyed. Destruction seems to be associated with the decelerating GAS spending 

significant time in the forbidden velocity region {Ps < \P\ < Pa-)- Figure [3] shows an unstable supersonic GAS, in 

which the soliton retains its integrity, settling to a stable subsonic GAS. 

Summary. — Electrostriction couples acoustic waves to light in a waveguide with a Bragg grating, and the resulting 

system supports "gap-acoustic" solitons, a generalization of optical gap solitons. Electrostriction damps the gap 
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soliton's (known) oscillatory instability, and gives rise to another instability, which occurs when the soliton's velocity 
exceeds the speed of sound. In contrast to the oscillatory instability, which destroys the gap soliton, the supersonic 
instability tends to leave the soliton intact, transforming it into a stable subsonic soliton. Thus, electrostriction adds 
a new level of complexity to the structure of optical gap solitons, and may facilitate the creation of ultra-slow optical 
pulses. Electrostriction and acoustic waves are present in virtually all materials, so an understanding of physically 
realistic optical gap solitons, especially slow ones, must entail a grasp of the effects of electrostriction. 
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